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Abstract 

The dynamics of particles interacting by key-lock binding of attached biomolecules are studied 
theoretically. Experimental realizations of such systems include colloids grafted with comple- 
mentary single-stranded DNA (ssDNA), and particles grafted with antibodies to cell-membrane 
proteins. Depending on the coverage of the functional groups, we predict two distinct regimes. In 
the low coverage localized regime, there is an exponential distribution of departure times. As the 
coverage is increased the system enters a diffusive regime resulting from the interplay of particle 
desorption and diffusion. This interplay leads to much longer bound state lifetimes, a phenomenon 
qualitatively similar to aging in glassy systems. The diffusion behavior is analogous to dispersive 
transport in disordered semiconductors: depending on the interaction parameters it may range 
from a finite renormalization of the diffusion coefficient to anomalous, subdiffusive behavior. We 
make connections to recent experiments and discuss the implications for future studies. 
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I. INTRODUCTION 



In this paper we present a theoretical study of desorption and diffusion of particles which 
interact through key-lock binding of attached biomolecules. It is becoming common practice 
to functionalize colloidal particles with single-stranded DNA (ssDNA) to achieve specific, 
controllable interactions ^1.31.8). Beyond the conceptual ,ntere.t as a .node, 
system to study glassiness[7|] and crystallization, there are a number of practical applications. 
Colloidal self-assembly may provide a fabrication technique for photonic band gap materials 
( 8|, 9|). One of the major experimental goals in this line of research is the self-assembly of 
colloidal crystals using DNA mediated interactions. The difficulty stems in part from the 
slow relaxation dynamics in these systems. The main goal of this paper is to understand 
how the collective character of key-lock binding influences the particle dynamics. In doing so 
we gain valuable insight into the relaxation dynamics, and propose a modified experimental 
setup whose fast relaxation should facilitate colloidal crystallization. 

Similar systems have also attracted substantial attention in other areas of nanoscience. 
In particular, by functionalizing nanoparticles with antibodies to a particular protein, 
the nanoparticles have potential applications as smart, cell-specific drug delivery vehicles 
( These nanodevices take advantage of the fact that certain cancerous cells over- 
express cell membrane proteins, for example the folate receptor. An improved understanding 
of desorption and diffusion on the cell membrane surface may have implications for optimiz- 
ing the design of these drug delivery vehicles. 

In what follows we present our results on the dynamics of particles which interact through 
reversible key-lock binding. The plan for the paper is the following. In section II we 
introduce the key-lock model and explain the origin of the two model parameters A and m. 
The parameter A determines the binding energy for the formation of a key-lock pair. The 
parameter m is the mean of the distribution for the number of key-lock bridges. Depending 
on m, which is related to the coverage of the functional groups (e.g. ssDNA), there are 
two distinct regimes. At low coverage there is an exponential distribution of departure 
times, but no true lateral diffusion. As the coverage increases, we enter a regime where the 
particle dynamics is a result of the interplay between desorption and diffusion. An estimate 
is provided for the value of m which determines the crossover from the localized to diffusive 
regime. In section III the localized regime is discussed in detail. In this regime the particle 
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FIG. 1: (Color online). Graphical depiction of particles interacting with a flat, 2D substrate by 
multiple key-lock binding. 

is attached to a finite cluster and remains localized near its original location until departing. 
We derive the partition function for the finite clusters, and calculate the departure time 
distribution. In section IV we determine the departure time distribution in the diffusive 
regime. We present an effective Arrhenius approximation for the hopping process and a 
Fourier transform method which greatly simplifies the calculation. In section V we discuss 
the random walk statistics for the particles' in-plane diffusion. A set of parametric equations 
is derived to relate the average diffusion time to the mean squared displacement. The 
lateral motion is analogous to dispersive transport in disordered semiconductors, ranging 
from standard diffusion with a renormalized diffusion coefficient to anomalous, subdiffusive 
behavior. In section VI we connect our results to recent experiments with DNA-grafted 
colloids. We then discuss the implications of the work for designing an experiment which 
facilitates faster colloidal crystallization. In section VII we conclude by summarizing our 
main results. 
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II. MODEL 



We now present the model, where a single particle interacts with a flat two-dimensional 
surface by multiple key lock binding (see Fig. [T]). At each location on the surface there 
are m key-lock bridges which may be open or closed, with a binding energy of e for each 
key-lock pair. Here we have neglected the variation in e. In the case of the DNA-colloidal 
system mentioned in the introduction, the model parameter e is related to the hybridization 
free energy of the DNA. The resulting m-bridge free energy plays the role of an effective 



local potential for the particle 



U{m) = -kBTmA (1) 
A = log(l + exp[e/fcBT]) (2) 

Generically, m is a Poisson distributed random number = fTf" exp(— m)/m! where fn 
denotes the mean of the distribution. The model parameter m is a collective property of 
the particle-surface system. For example, consider the case of dendrimers functionalized 
with folic acid, which can be utilized for targeted, cell speciflc chemotherapy. The folic 
acid on the dendrimer branch ends form key-lock bridges with folate receptors in the cell- 
membrane. In this case fn will depend on the distribution of keys (folic acids) on the 
dendrimer, and the surface coverage of locks (folate receptors) in the cell membrane. 

At each location, the particle is attached to the surface by m bridges. To detach from the 
surface the particle must break all its connections, in which case it departs and diffuses away 
into solution. Alternatively the particle can hop a distance a to a new location characterized 
by a new value of the bridge number m. By introducing the correlation length a, we have 
coarse-grained the particle motion by the distance after which the new value of the bridge 
number becomes statistically independent of the value at the previous location. In the 
localized regime the particle remains close to its original location until departing. In the 
diffusive regime the particle is able to fully explore the surface through a random walk by 
multiple breaking and reforming of bridges. 

Naively one might expect the crossover between the two regimes to occur at the percola- 
tion threshold, where one flrst encounters an inflnitely connected cluster of sites with m > 0. 
However, the crossover from the localized to diffusive regime occurs at smaller m than pre- 
dicted by percolation theory. If pc = 1/2 denotes the critical probability for site percolation 
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FIG. 2: (Color online). The crossover from the localized to the diffusive regime below the per- 
colation threshold. Estimates based on the characteristic cluster size (Eq. [3l dashed line) and 
confinement of the random walk (Eq. HI solid line) give similar crossovers. For large A the 
crossover condition is m = — log (l — 4) . 



on the triangular lattice, the percolation transition occurs at m = log (2). There are two 
alternative estimates for the crossover from the localized to the diffusive regime. The first 
is to compare the average number of steps n = exp (Am) the particle takes before departing 
(see section IV) to the characteristic cluster size Sc = l/\og{l/\p) below the percolation 



threshold. Here A = 5.19 is a numerical constant for the triangular lattice [13|], and in 
the percolation language p = 1 — exp(— m) is the occupancy probability. The crossover 
condition n = Sc can be expressed as a function of m. 



A = --log [-log{A(l-e-)}] 



(3) 



Alternatively, in the localized regime the particles' random walk is confined by the charac- 
teristic cluster size. Below percolation the radius of gyration of the cluster is Rs ~ with 
p = 0.641 in two dimensions. Comparing the radius of gyration of the cluster to the radius 
of gyration for the particles' random walk, the crossover occurs at n = (sc)^''. 

2p 



A 



m 



log [-log{A(l-e-'")}] 



(4) 



Since 2p differs from 1 by less than 30%, both conditions give similar crossovers (see Fig. 
[2]). The saturation at m = — log (l — ^) occurs for very large A, as a result for binding 
energies of a few ksT per bridge the crossover occurs at m ~ 0.1. 



III. LOCALIZED REGIME 



In the percolation language, when the occupancy probability p = 1 — exp(— m) is small, 
particles are localized on finite clusters. In this localized regime particles are able to fully 
explore the cluster to which they are attached before departing. This thermalization of 
particles with finite clusters permits an equilibrium calculation of the cluster free energy F = 
—ksTlogl^Z) . The departure rate is given by the Arrhenius relation K = exp (F/kBT). 
Here tq is a characteristic timescale for bridge formation. The probability that the particle 
departs between t and t + dt is determined from the departure time distribution ^{t)dt ~ 
Kexp[-Kt]dt. 

To begin we calculate the partition function for the finite clusters. The cluster is defined 
as s connected sites on the lattice, all of which are characterized by < m < m* bridges. 
For Poisson distributed bridge numbers the partition function for the finite cluster is: 

s m*~l 

Z{m*, s) = y^ Pm, exp(Ami) = — -(exp(me^)Q(me'^, m*) - 1) (5) 

^ exp(m) — 1 

Because by definition the cluster does not contain sites with m = bridges we have 
renormalized the probability distribution Pm = Pm/{^ — exp(— m)) so that J2m=i^m = 
1. Here Q{x,m*) = T{x,m*)/T{m*) = exp{—x)J2T=o^^'^/^^- regularized upper 

incomplete F function. In the language of the statistics of extreme events, m* — 1 is the 
maximum "expected" value of m in a sample of s independent realizations [l^. The point 
is that on finite clusters we should not expect to achieve arbitrarily large values of the 
bridge number. Hence when averaging the partition function to obtain the cluster free 
energy one should only average over sites with m < m*. The distribution function for 
m* is obtained by noting that the probability that all s values of m are less than m* is 
{Ylm=i^ Pm j = (^ ^^"^^^xpim)'-! ^~^ ) ■ differentiating this quantity with respect to m* we 
obtain the distribution function for the maximum expected value of m. 

(6) 

\ exp(m) — 1 / 
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FIG. 3: (Color online). Departure time distribution function versus time in the localized regime 
with m = 0.05. The results of our calculation (solid lines) are compared to single exponential 



relaxation with departure rate K 



exp(-A) 



(dotted lines). 



The cluster size distribution below the percolation threshold is exponential [1^] with char- 
acteristic cluster size Sc = l/log(l/Ap). 



1-Xp f s 
Ps{p> = — 7 — exp 

A \ Sr 



(7) 



The summation over s can be performed analytically, which allows the result to be expressed 
as a single summation over m*. 



CO oo 



(8) 



m*=2 s=l 



A(l - A(l - exp(-m))) 



exp(m) — 1 



Pm*(exp(me^)Q(me^,m*) — 1) 



m*=2 

y{m*) = X{Q{m,m*) — exp(— m)) 



1 + y{m*) 
(1 — y{in*)Y 



(9) 



For fixed m, as in the plot (see Fig. [3]), changing A is directly related to a change in 
the average binding free energy. Increasing A leads to a reduction in the rate of particle 
departure. 
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IV. DIFFUSIVE REGIME 



The departure time distribution changes significantly in the diffusive regime. In this 
regime the particle can explore the surface to find a more favorable connection site, which 
leads to a longer lifetime for the bound state. This phenomenon is qualitatively similar 
to aging in glassy systems. In these systems one finds that the response to an external 
field is time dependent [16]. In the magnetic analogy this leads to a time dependence of 
the magnetization. Below the glass temperature, the longer one waits before applying the 
external magnetic field, the more time the system has to settle into deep energy wells, and 
the smaller the response. In our case, the diffusive exploration of the particle allows it to 
find a deeper energy well, which leads to an increase in the bound state lifetime. 

As a result, the departure time distribution must now reflect not only desorption, but 
also hopping to adjacent sites. The hopping rate between neighboring sites i and j is given 
by an Arrhenius law Ki^j = :^exp[— A(mj — mj)6{mi — rrij)], with 6{x) the Heaviside step 
function. In a lattice model with coordination number z the dwell time at a site with m 
bridges is calculated by averaging over the hopping rates to the nearest neighbors (see Fig. 

m). 

/1\ °° °° ^ ^ 1 

^^ = ( iv^. ) =^E---E^----^-v^^^ — (10) 

Fortunately, this ensemble averaging procedure can be accurately approximated by an ef- 
fective Arrhenius relation: 

= roexp[A(m-m)] (11) 

The validity of the approximation is most important for sites with m >m bridges, since the 
diffusive exploration allows the particle to quickly cascade into these deep energy wells. 

The effective Arrhenius relation greatly simplifies the calculation, since so long as Am is 
sufficiently large, the probability of the particle still being attached to the surface after an 
n step random walk is (1 — KmTm)"'~^ — ^ exp(— Am)]"""'^. Here = ^ exp{—Am) is 
the departure rate from a site with m bridges. Interestingly, in this approximation scheme 
the attachment probability is independent of the particular bridge numbers {mi,...,m„} 
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m 

FIG. 4: (Color online). Comparison of the ensemble averaged dwell time in a lattice model (Eq. 
[To]) to the effective Arrhenius approximation (Eq. [TT|) . In the plot the product Am = 4 is held 
constant. 

realized during the walk. Thus, the probability of departure /„ after exactly n steps is: 

/„ = [exp(7) - 1] exp(-7n) (12) 
7 = — log[l — exp(— Am)] (13) 

The average number of steps for the random walk is "^^^=1 ^fn = exp(Am). To calculate the 
departure time distribution $(t) we use /„ to average over the departure time distribution 
for walks with a given n, (j)n{t)- 

oo 
n=l 

M^)-u[±^.jy,{=^))s{.-p)l (15) 

Here Sm{t) is the survival probability at time t for a site with m bridges, used to determine 
the probability of departure between t and t + dt. If there was only one hopping pathway 
with rate k, we would have — ^ = Kexp(— Kt). The generalization accounts for the fact 
that the particle can hop to any of its z neighbours, and the probability of departure is not 

9 



simply exponential. 



(16) 



,a=l 



It is convenient to Fourier transform so that one can sum the resulting geometric series 
for 



/oo 
(f)n{t) exp[-iut]dt = Xiu)"" 
■oo 



X{u;) = ^ J] • • • 5^ P^, ■ ■ -P, 



m=l mi=l mz=l 



1=1 



i=l 



+ iuj 



$(a;) = [exp(7) - 1] [exp( -7)^(0;)]" = [exp(7) - 1] 



Xiuo) 



(17) 
(18) 
(19) 



„=i exp(7) - X{uj) 

To facilitate a simpler calculation, we employ a coarse-graining procedure to dispense 
with the tensor indices {m,mi, ...,777.2} in the definition of X{ijj). In the summation there 
are many terms for which the value of ^ Hm^mi equal, but with different weight fac- 

^ ^ ^ i=l * 

tors PmPmi ■ ■ 'Pruz- To eliminate this degeneracy we introduce a smooth function /(k) 
normalized according to J f{K)dK — 1. 



X(uj) ~ / fU)^^dK 
J K + IUJ 



(20) 



The inverse Fourier transform is performed using the residue theorem to obtain the final 
result. The contour integral is closed in the upper half plane, with all the poles on the 
imaginary axis ai lu = iz. 

exp[iujt]X {uj) 



27r 



27r 



r=l 



exp(7) -X{uj) 



= [exp(7) - l]i ^ res, 

r=l 

00 

= [exp(7) - l]i 



r=l 



exp(7) - {Xiuj^) + {uj- ujr) + • • • } 

exp[iuJrt]X (cUr) 



(dX) 



exp(7)[exp(7) " 1] ^ ^T^(^^''^^ 



Y{zr) = / /(«:) 



K 



{k — ZrY 



(21) 



(22) 
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FIG. 5: (Color online). Departure time distribution function versus time as determined by Eq. 
[2T] in the diffusive regime. In the plot the average binding energy is held constant at SksT. 
The theoretically determined departure time distribution can be compared to an experiment with 
DNA-grafted colloids (see Section VI) which observed power law behavior with exponent —1.5. 

Here labels the roots of the equation 

exp(7) - X{iz) = (23) 

The benefit of the coarse-graining is now more transparent, as the residues are all labeled 
by a single index r as opposed to the tensor indices {m, mi, m^}. 

In Fig. [5] the departure time distribution is plotted in the diffusive regime. The optimal 
regime for fast particle departure is to have a large number (m ~ 10) of weakly bound key- 
lock bridges. In this scenario the departure time distribution is accurately approximated 
as a single exponential, ^(t) = Kmexp{—Kmt). 

V. DIFFUSION 

We now turn to discuss the statistics for the in-plane diffusion of the particle. We 
first note that the in-plane trajectory of the particle subjected to a delta-correlated random 
potential remains statistically equivalent to an unbiased random walk. As a result, the 
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FIG. 6: (Color online). Root mean squared displacement vs. time with Am = 4. The curves are 
calculated from the parametric equations [261 123 



mean squared displacement for an n step random walk remains (r^) 



na 



As the particle 



explores the landscape it cascades into deeper energy wells, the hopping time increases, and 
the diffusion gets slower. In the limit n oo the average hopping time can be determined 
from the equilibrium canonical distribution. For Poisson distributed bridge numbers m, this 
corresponds to a finite renormalization of the diffusion coefficient D* with Dq = /Atq. 



(t) = n {Tm) = utq exp (—Am) ^ Pm exp (Am) = utq 

m=l 



exp (me^) 



1 d (r^) 
D* = = Do 



exp (Am) [exp(m) — 1] 
exp (Am) [exp(m) — 1] 



(24) 



(25) 



A d{t) " exp (me^) - 1 
This "ergodic" behavior is only achieved after a very long time. Generally, an n step 
random walk cannot visit sites with arbitrarily large m. In this transient regime one should 
only average over sites with m < m*. In the language of the statistics of extreme events, m* — 
1 is the maximum "expected" value of m in a sample of n independent realizations [14]. Even 
with this complication, the average diffusion time (t) and the mean squared displacement 
(r^) can both be expressed in terms of m*, which defines their relationship in parametric 
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1 dlr'^) 

FIG. 7: (Color online). The dimensionless diffusion coefficient D = 475^-^^ plotted against 
time. 



form. 



P{m, m*) 



(26) 



(27) 



D* \ 1 — exp(— me^^ 
Here P{x,m*) = 7(x, m*)/r(m*) = exp(— x) x'^/fc! is the regularized lower in- 

complete r function. In the limit m* ^ 00 we recover the renormalized diffusion relation 
(t) = (r^) /D*, although this occurs at very long, often unrealistic times. In the transient 
regime we expect anomalous, subdiffusive behavior. As indicated in Fig. El this subdiffusive 
behavior is typical for strong enough key- lock interactions. Figure [7| is a plot of the dimen- 
sionless diffusion coefficient versus time. The anomalous diffusion may be well described 
by a power law (r^(t)) ~ (t)" with a nonuniversal exponent a < 1. 

There is an analogy between our results and dispersive transport in amorphous materials. 
In these systems a length dependence of the effective mobility[l7| can be interpreted within 
the context of the statistics of extreme events The time required for charge carriers 

to travel through the material depends on the dwell times spent at all of the trapping 
centers. Since this transit time will be dominated by the dwell time of the deepest trapping 
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center, one would like to know how the thickness of the material effects the distribution for 
the largest trapping depth. The analogy to our result is made by replacing the material 
thickness with the number of steps in the random walk n, and replacing the distribution for 
the largest trap depth with the distribution for m*, since the bridge number is related to 
energy by Eq. [H 



VI. EXPERIMENT 



We now would like to make a connection between our results and recent experiments with 
DNA-grafted colloids. The departure time distribution can be compared to an experiment 
which determined the time-varying separation of two DNA-grafted colloids in an optical 
trap l|. In the experimental setup, two particles are bound by DNA bridges, and after 
breaking all connections diffuse to the width of the optical trap. Because the length of 
the DNA chains grafted on the particle is much shorter than the particle radius, surface 
curvature effects can be neglected. The interaction resembles that of a particle interacting 
with a patch on a 2D substrate. Experimentally the tail of the departure time distribution 
was observed to be a power law ~ Qualitatively similar behavior is predicted by 

the theory with m ~ 1 and average binding free energy of several ksT (see m = 1 curve in 
Fig. ED. 

In addition, our work provides insight into the slow crystallization dynamics of key-lock 

binding particles (see Fig. [8]). In [2j|, l/im diameter particles grafted with ssDNA formed 

reversible, disordered aggregates. The average number of key-lock bridges between particles 

n 

was m ~ 2. The authors of [1] observed random hexagonal close packed crystals by further 
reducing the surface density of DNA strands on the particles. The crystallization process 
requires that particles rapidly detach and reattach at the desired lattice location. In the 
localized regime particle desorption is the relevant process. 

In the diffusive regime surface diffusion also plays a role in the rearrangement of particles 
into the desired crystalline structure. To determine which process is more important for 
particle rearrangement, we can compare the departure time with the time required for a 
particle to diffusively explore the surface of a particle to which it is bound. The time T^ep 
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FIG. 8: (Color online). Graphical depiction of key- lock binding between nanoparticles func- 
tionalized with complementary ssDNA. The resulting structures can be disordered, fractal-like 
aggregates, or crystalline. 




FIG. 9: (Color online). Plot of the characteristic times T^^gp and T^if versus m at constant binding 
energy (const = 4 in Eq. [29]) . 
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required for 90% of the particles to depart is: 



0.1 = / <^{t)dt (28) 

To estimate the time required for diffusive rearrangement Tdif we use the parametric 
equations In particles of radius R = .5/im were grafted with DNA chains of 

length / ~ 20nm. Assuming the correlation length a ~ / we have ~ (t^) — -'-O^- 
Figure [9] shows a comparison of Tdif and Tdep at constant binding free energy. 

/\Tn 

, . + log(l — exp(— m)) = const (29) 

1 — exp(— m) 

This expression for the binding energy takes into account the entropy reduction associated 
with the non-ergodic degrees of freedom. For a detailed discussion of this topic see reference 



12l |. Since T^ij > T^ep, colloidal desorption and reattachment is the dominant mechanism 



by which particles rearrange. 

As the figure indicates, the optimal regime of fast departure is to have a large number 
(m > 10) of weakly bound key-lock bridges. We predict a localized regime below the 
crossover where particle departure is relatively fast. Just beyond the crossover there is a 
relative maximum in T^ep before it decreases at large m. The increase in departure time 
at the onset of diffusive behavior is indicative of a regime where the system ages. In this 
regime the interplay of diffusion an desorption leads to longer bound state lifetimes, and an 
increase in the departure time. 

We now turn to the question of designing a future experiment which will facilitate fast 
particle departure and hence colloidal crystallization. The essential goal is to increase the 
average number of key-lock bridges between particle pairs. Increasing the surface density 
of DNA strands alone results in the formation of a brush, which decreases the effective 
cross section for the interaction between complementary DNA. Instead we propose the 
introduction of long, flexible DNA linkers between particles with a high coverage of short 
ssDNA (see Fig. ITOl) . This system has the potential to realize more key-lock bridges between 
particle pairs as compared to previous experiments, and therefore substantially reduce the 
time required for crystallization. 
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FIG. 10: (Color online). Graphical depiction comparing recent experiments with DNA-grafted 
colloids to a future implementation with long flexible linkers. Increasing the number of key-lock 
bridges between particle pairs potentially decreases the time required for crystallization. 

VII. CONCLUSION 

This paper studied the dynamics of particles which interact through the reversible forma- 
tion of multiple key-lock bridges. Well before the percolation threshold is reached there is a 
crossover from a localized regime to a diffusive regime. In the localized regime the particles 
remain close to their original attachment site until departing. In this regime particles are 
attached to finite clusters, and the system exhibits an exponential distribution of departure 
times. Once the radius of gyration of the cluster exceeds the characteristic radius for the 
particles' random walk, the finite clusters behave effectively as infinite clusters. Diffusion 
allows the particles to cascade into deeper energy wells, which leads to a decrease in hopping 
rate. The diffusion slows and the bound state lifetime increases, a phenomenon qualitatively 
similar to aging in glassy systems. In the diffusive regime we discussed the statistics for the 
particles' in-plane diffusion. Weak key-lock interactions give rise to a finite renormalization 
of the diffusion coefficient. However, as the strength of the interaction increases (larger 
A), the system exhibits anomalous, subdiffusive behavior. This situation is analogous to 
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dispersive transport in disordered semiconductors. We then made the connection between 
our calculation of the departure time distribution and recent experiments with DNA-coated 
colloids. The findings indicate that the optimal regime for colloidal crystallization is to 

have a large number of weakly bound key-lock bridges. A modified experimental setup was 
proposed which has the potential to realize this regime of fast particle departure. 
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